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ABSTRACT 

We constrain the possible time variation of the Higgs vacuum expectation 
value (v) by recent results on the primordial 4 He abundance (Yp). For that, we 
improve the analytic models of the key-processes in our previous analytic calcu- 
lation of the primordial 4 He abundance. Furthermore, the latest results on the 
neutron decay, the baryon to photon ratio based on 5-year WMAP observations 
and a new dependence of the deuteron binding energy on v are incorporated. 

Finally, we approximate the weak freeze-out, the cross section of photo- 
disintegration of the deuteron, the mean lifetime of the free neutron, the mass 
difference of neutron and proton, the Fermi coupling constant, the mass of the 
electron and the binding energy of the deuteron by terms of v, to constrain 
its possible time variation by recent results on the primordial 4 He abundance: 
|^| < 1.5 ■ 10" 4 . 

Subject headings: cosmology: theory — cosmology: cosmological parameters - 
cosmology: early universe 



1. Introduction 



The standard model (lGriffithslll987l ) is a remarkably successful description of fundamen- 
tal particle interactions. The theory contains parameters - such as particle masses - whose 
origins are still unknown and which cannot be predicted, but whose values are constrained 
through their interactions with the conjectured Higgs field. The Higgs field is assumed to 
have a non-zero value in the ground state of the universe - called its vacuum expectation 
value v - and elementary particles that interact with the Higgs field obtain a mass propor- 
tional to this fundamental constant of nature. 
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Although the questi on whether the fundamental constants are in fact constant, has a long 
history of study (see lUzanl (120031 ) for a review), com paratively less interest has been directed 



towards the consequences of a possible variation of v (jGaBner fc Le sch 200 81; iDent et al 



2002; 



Kuiat & Scherreriboool : IScherrer & Spergell Il993l : bixit & Sher1ll988h . 



2008; 



Chamoun et~al1l2007l:Eandau et al.lbood: IlI fc Chull20od:lYoo j ScherreJ2003l : llchikawa fc Kawasaki 



A macroscopic probe to determine the allowed variation ran ge is given by th e network of 
nuclear interactions during the Big-Bang-Nucleosynthesis (see lYao et al.l (120061 ) for a review 
of the Standard Big-Bang- Nucleosynthesis Model SBBN), with its final primordial abun- 
dance of 4 He. The relevant key-parameters are the freeze-out concentration of neutrons and 
protons, the so called deuterium bottleneck (the effective start of the primordial nucleosyn- 
thesis) and the neutron decay Their dependency on v and the final impact on the resulting 
primordial abundance of 4 He can be understood more clearly by an indepth approach. 



Here we present a revised calculation of the primordial 4 He abundance (jGaBner fc Lesch 



20081 ). where the analytic models of all key-processes have been improved. The opening of 
the deuterium bottleneck and the weak freeze-out are determined more accurately and a 
new dependence of the deuteron binding energy on v is incorporated, based on different 
nucleon-nucleon-potential- models. 

The analytic approach enables us to take important issues into consideration, that have been 
ignored by previous authors, as there are the f-dependence of the relevant cross sections of 
deuteron production and its photo-disintegration. Furthermore, we take a non-equilibrium 
Ansatz for the freeze-out concentration of neutrons and protons and incorporate the latest 
results on the neutron decay and the baryon to photon ratio. 

Finally, we approximate the weak freeze-out, the cross section of photo-disintegration 
of the deuteron, the mean lifetime of the free neutron, the mass difference of neutron and 
proton, the Fermi coupling constant, the mass of the electron and the binding energy of 
the deuteron by terms of v, to constrain its poss i ble time variation by recent results on the 
primordial 4 He abundance (IPeimbert et al.l 120071 ; llzotov et al.l 120071 ). 



We briefly note, that constraints on the spacial variation of v require a measurement 
of helium abundance anisotropy or inhomogeneity versus the position in the sky and an 
inhomogeneous theoretical BBN model. The homogeneous formalism used throughout the 
paper thus assumes a spacial invariance of the Higgs vacuum expectation value. 
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2. Calculations 

All relevant processes of SBBN took place at a very early epoch, when the energy density 
was dominated by radiation, leading to a time-temperature relation for a flat universe: 



1 r i m 

327r 3 k 4 Gg*T 2 [S1, U 

where c is the velocity of light, k the Boltzmann constant, G denotes the gravitational 
constant and h is the Planck constant divided by 2tt. g* counts the total number of effectively 
massless (mc 2 <C kT) degrees of freedom, given by g* = (g^ + I*?/), in which represents 
the bosonic and gf the fermionic contributions at the relevant temperature. 
At very high temperatures (T ^> 10 10 K), the neutrons and protons are kept in thermal and 
chemical equilibrium by the weak interactions 

n + e + p + & e , 
n + v e p + e~ and 
n ^ p + e" + u e , 

until the temperature drops to a certain level, at which the inverse reactions become ineffi- 
cient. This so called "freeze-out" -temperature Tf and time tf denote the start of the effective 
neutron beta decay. 

Assuming chemical and thermal equilibrium, the rate of neutron to proton concentration 
at freeze-out is commonly calculated, assuming chemical and thermal equilibrium: 

-(?))= e~*, (2) 

Tip 

where Q denotes the energy difference of neutron and proton rest masses. However, the 
deviation from equilibrium at freeze-out is already significant. Hence, we have to der i ve no n- 



equilibrium concentrations, where we follow the example calculations of iMukhanovl (J200J): 
The 4-fermion-interaction a + b —>■ c + d can be calculated using the Fermi theory, where the 
differential cross section is given by 



dogb _ 1 |Mf / (go • Pd) 2 - m 2 c m\ c 8 
dVt (8tt) 2 (p a +p b ) 2 \j {pa-Pb) 2 -m 2 a ml c 8 ' 

{Pa ' Pb) an d (p c ' Pd) denote the scalar products of the 4-momenta and the matrix element is 
given by: 

\M\ 2 = 16(1 + 3g 2 A ) G 2 F (p a -Pb)(pc-Pd), 
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where G F ~ 1.166 371 ■ 10 



-ii 



and g A = 1.2739 (lAbele et al. 



McV 



(jYao et al.ll2006l ) denotes the Fermi coupling constant 



20021 ) the axial vector coupling constant, respectively. 



First, we consider the reaction n + v e — > p + e at the relevant temperatures around few 
MeV and below, where the nucleons are nonrelativistic: 



On + P. 
(Pn ■ Pv) 
(Pp ■ Pe) 



2 2 4 

~ m n c 



m n c e v 
m p c 2 e e 



yj (Pp ■ Pe 



) 2 — m^mgC 8 



m p c e e v e 



where m p , m n and m e denote the mass of the proton, neutron and electron, respectively, v e 
is the velocity of the electron, e u is the energy of the incoming neutrino and e e ~ e v + Q is 
the energy of the outgoing electron. Substituting all terms into eq. we obtain: 



da r . 



(Svr) 2 
1 



16(1 + 3^ A )G F 



m n a/ im n t v f - (m n m u c 2 y 



,1 + %9a) G 2 F — | e 2 — , 
m„ c 



(4) 
(5) 



where we neglect the neutrino mass m v . Integration leads to 



Or. 



l + 3 5 



7T 



A n 2 m p 2 V e 

m 2 n c 



(6) 



Next, we have to consider that at temperatures kT > 2 m e c 2 the possible states for the 
electron are partially occupied by electron-positron-pairs. According to the Pauli exclusion 
principle, this reduces the appropriate cross section to 



a 



nu __£e_ 
1 + e kT 



(7) 



This enables us to calculate AN n , the reduction of neutrons within a time interval At in a 
given volume, containing N n neutrons: 



AN„ 



< u n ev v„Ag eu N n At. 



(8) 



where 



n r 



1 

1 + e& 



(9) 
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denotes the neutrino occupation number (v u and T v are the velocity and the temperature of 
the neutrinos) and 



&9e 



2tt 2 



<E„ + A<!„ 



\P? d\p\ 



Z7T 



the phase volume element. Introducing the relative concentration of the neutrons 



A. 



n.,. 



n n + n p 



(10) 



(11) 



and assuming baryon conservation, we obtain the rate of change of the neutron concentration 
due to the ni/-process: 

'dX„ 



dt 



^■nu-^-ni 



(12) 



where X nu denotes the decay rate. Substituting the cross section (JTj) into eq. (JHJ) we obtain: 



L+jgi C 2 m P 1(T s 



where 



3{T V 



el Wl 



m e c 2 \ 2 \ 



1 + e kT x + ekTv 



— de l 



(13) 



(14) 



Below the temperature kT ~ 2 m e c 2 , the Pauli exclusion principle, represented by the 
term (1 + e~fcr), looses importance and numerically we notice a deviation of 1 % only, when 
we set this term to 1. Expanding the square root (m e c 2 /e e <C 1) keeping only first two terms 
and introducing the integration variable x = we derive: 



3(T„ 



(kT v 



X 



Q 



1 / m e c J 



2 \ kT v 



1 + e x 



dx 



(15) 



kT, 

Q 



45 C(5) (kT, 



Q 



7n 4 (kT 



60 V Q 



Q 



5 45 C(5) fkTA 6 fkT, 

2 V Q J V Q 



+ 0.25 



2 4 



303) / 

2 V 2Q 2 



(16) 



(17) 



where ( is the Riemann zeta function. In the last step, we completed the square approxi- 
mately. Finally we convert X nu from MeV to - and derive 



A ~ 



l + 3flj r2 fkTA 3 (m 

1.75 ■ 10~ 21 fn 2 \ Q J V Q 



+ 0.25 
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Similarly, we find the decay rate of the reaction n + e + — > p + v (we interchange with e e 
and m e with m v = 0): 

Ane = 7T~r^ — f / e l 1 , 6 e e <^ e (19) 

2tt 3 m2 J mec2 1 + efcr 

Assuming T„ = T, the rates of the inverse reactions are related to the rate of the direct 
reactions as 

Ap e = e'^Kv (20) 

\v = e~*r\ ne . (21) 

Hence, we can write the following balance equation for X n : 

dX 

- ~{^nu + A ne )X n + (Ap e + \ pv )\l — X n ) 



dt 

= -(\ nv + \ne)(l+e-£)(X n -X^) (22) 

with the equilibrium neutron concentration 

K q = —^o- (23) 

1 + ekT 

To solve this linear differential equation ( 1221) . we take the initial condition X n (t = 0) = X^ q 
and obtain: 



X n (t) = X?(t) - J exp f - y_ (X nv (y) + X ne (y))(l + e~^)dy )X?(t)dt, (24) 

where dot denotes the derivative with respect to time. Using the auxiliary function F(t) 

F(t)= / (\ w (t) + \ ne (t))(l + e-&)dt, (25) 
Jo 

we express the integral in ( 1241) in the form 

[ l e -F(t)+F(i) di (26) 

JO 

and expand in the small parameter (t—i), since the integral is dominated by the contribution 
of t ~ t if F(t) is a quickly growing function of t: 

Xnit) = X«(t) - jf [X«{t) + K\t)(t-t) + ...] e~*M (l + l -F(t)(t - tf + dt(27) 
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We integrate term by term using 

ft 



-A(t-t) 



(t - t) n dt 



(28) 



o 



where we neglect exponentially small terms of order e , deriving: 



Xn(t) = X?(t) 1 



(A^ + A ne )(l + e-*)^(t) 



(29) 



Later, when the temperature has dropped significantly, X^ goes to zero and the integral 
in eq. (Elj) approches the finite limit. As a result, the neutron concentration freezes-out 
at X n (t — > oo). Effectively, this freeze-out occurs, when the deviation from equilibrium 
becomes significant, hence when 



X eq 

n 
An 



{Kv + A ne )(l + e w). 



(30) 



Assuming this happens before e ± -annihilation and after kT has dropped below Q, we set 
A™, + A„ e ~ 2A nj , and neglect the term exp(— Q/kT). Substituting all terms and taking the 
time-temperatur- relation ([1]) into account, we finally obtain an equation for the freeze-out- 
temperature Tf\ 

2 



G9 * {Tf) 3.7-lQ- 35 = (l + 3g 2 A )G 2 F Q 3 f^f-Y ( ^ + 0.25 ) . (31) 



Q 



Q 



The quadratic term ^S- leads to 



6 Q 



T f ~ 1.16 • 10° ^ 



1 



\ 



3.7 • 10" 35 



1 

64 + V (1 + 3^) G 2 F Q ? - 



Gg*{T f ) 
h 3 c 5 



1/4 



[K]. (32) 



At Tf the effectively massless species in the cosmic plasma are neutrinos (left-handed only), 
antineutrinos (right-handed only), electrons, positrons and photons. For the case of three 
neutrino families (N„ = 3), we obtain g*(T f ) = (2 + f (4 + 2N U ) = 10.75. 
To calculate the relevant neutron concentration at Tf, we go back to eq. ([241 . Since X^ q — > 
as T — > we have to calculate the integral term in eq. (1241 in the limit t — > oo. The main 
contribution to the integral comes at temperature above the restmass of the electron, where 
again A n „ + X ne ~ 2 X nu (X nu given by eq. (JIB"]) ). Furthermore, we use eq. ([T]) to change the 
integration variable from to dT: 



Xn(Tf) 



Q exp 



-2.68 ■ 10 



31 



Os.(r/) 



2 T 2 (1 + cosh2) 



dT, (33) 
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with Q and T in units of MeV. For later purposes, we finally state the neutron to proton 
ratio at freeze-out: 

= r (34) 

In comparison, the equilibrium rate eq. ([2]) is 14 % higher, thus as mentioned before, the 
deviation is significant, which justifies the effort. 

From now on th e loss of free neutrons via n — > p + e~ + u e , with a mean lifetime 



( ISerebrov et al.l 120051 ) r n = 878.5 s, can no longer be compensated. Thus, whereas the 

t-t i 

neutron density decreases as n n (t) = n n (tf) ■ e~~^~ , the proton density increases as n p (t) = 
n p (tf) + (n n (tf) — n n (t)) and we obtain 

n n X n (T f ) 

At the relevant densities in the early universe, fusion reactions can only proceed ef- 
ficiently through sequences of two-body-collisions and the starting product of these colli- 
sions is the weakly bound deuteron (B^ ~ 2.225 MeV), which is highly affected by photo- 
disintegration. The start of nucleosynthesis, tjsr, is therefore usually referred to as the "deu- 
terium bottleneck" . 

Once the deuteron production dominates the photo-disintegration and the expansion 
of the universe, our calculation in some sense "produces only deuteron", disregarding that 
deuteron is also destroyed by the fusion of light elements. In fact, we do not consider 
the detailed fusion reactions with their intermediate products that finally lead to 4 He. We 
are interested in the point t^, from that onwards we can assume neutron conservation, 
because enough neutrons have reached stable states inside light nuclei (no matter whether 
inside deuterons or further fusion products of deuteron). Hence, we obtain assuming two 
constraints: First, the deuteron production must dominate the photo-disintegration and the 
expansion of the universe. Second, to justify neutron conservation, the deuteron density 
must exceed the density of free neutrons. In fact, it turns out, that we can reproduce the 
numerical as well as the observational results for Yp, assuming neutron conservation when 
52 — 66 % of the neutrons have reached stable states inside light nuclei. 
The interval between tf and is a substantial fraction of the neutron lifetime and therefore 
plays an essential role for the outcome of the primordial helium production. Thus, we have to 
calculate the rates of deuteron production T^ np ^. dj ), deuteron photo-disintegration T^ d ^ np ) 
and the rate of deuteron density reduction by the expansion of the universe T expansion , to 
determine fx respectively Tn, when 

^ > {np—*d , y) ^ I '(-yd- *np) ^expansion- (36) 
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The rates for production and photo-disintegration of deuteron are given by the product of 
the relevant number densities, velocities and cross sections, whereas the expansion rate of a 
radiation-dominated, flat universe is given by ^ with t from eq. (jTJ, leading to: 



n„ 



77 n 7 

I 1 "a 



8kT 



rid 



V(np^d~f) > Ud C a^d^np) + 2~7 ; 



(37) 



where (T^d-mp) denotes the cross section of deuteron photo-disintegration, <7( np _»d 7 ) the cross 
section of deuteron production , tun is the nucleon m ass, 77 ~ 6.226 • 10~ 10 is the baryon to 
photon ratio based on WMAP (IKomatsu et al.ll2008l ) and lSteigmanl ( 120061 ). rid and n 7 denote 
the number densities of deuterons and photons, respectively, and n* the number density of 
photons which supply enough energy to disintegrate the deuteron and do not loose this 
energy in much more likely Compton scattering on electrons. 
The number density of photons at a certain temperature T is given by 



■in 



(he) 



ekT — 1 



"kT\ 3 
-i\E- = 16vr C(3) ( — J 



(3* 



where ( is the Riemann zeta function. The number density of these photons supplying a 
minimum energy E 1 > Bj 3> kT is 



n h>B d ) 



8vr 



Et t e~^dE 1 



57T 



kTX 
~hc~) 



Bd 
kT 



+ 1 



1 



e kT ? 



but most of them will loose energy in Compton scattering on electrons, leading to 

ild ^( / yd—>np) 



n 7 = n^ >Bd) 



rip (X("fe—>e"/) 



(39) 



(40) 



where cj( 7e ^ e7 ) denotes the Klein-Nishina cross section (IRybicki fc Lightman!ll979l ) for Comp- 
ton scattering on electrons: 



^(-ye— >e7) 

where 



7T \€om e c' 



1+/3 /2(l + /3) In (1 + 2/3) \ In (1+2/5) 1 + 3/3 



1 + 2(3 



ft 



+ 



2(3 



(i + 2py 



(3 



m e c 2 



(41) 



(42) 



and the mean incident photon energy (E^) is given by 
1 



•>7T 



n ii>B d ) (he) 



Ba 



E^-^dE^ 



kT 



\kT) 



B d + kT. (43) 
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The interaction cross sec tion of deuteron photo-disintegration can be well approximated 
by (IBethe &: Longmirdll950l ): 



<Jbl = El + Ml (44) 
where we express the electric dipole contribution 



El = I ^ft^-flJ- m 
3 ce m N (1 - r -±^/m N B d ) 



and the magnetic dipole contribution 



e 2 h(flp _ j Bd , Bd y (l - yfr^Eft + a.(r, + r t )g*ft - « s (r 5 - r f ) 



2 

m N (E 1 -B d ) s " 



poo 

-i2 



<X( 7 <^ np ) = 7T~T3 / E 7 e fcT <W# 7 ) (47) 

n (l>B d ) \ nc ) 



We go back to eq. (13T)|) and divide by rV 7d _,. n p), in order to receive two terms, which we 
analyse separately: 



r (np^dy) _ -. ^expansion 



1 1 expansion (48) 



(■yd^np) 1 (-yd^np) 



We start with 



l^d— »np) 



0"(7d— mp) 



(l + S) [(§ + l) 2 + l 

where (T( np _ > rf 7 ) is related to o r ( 7d __ >np ) by the detailed balance 

3(£ 7 > 2 

a^ d ^ np ) 2m N c 2 {(E 1 ) - B d ) ' 



(49) 



(50) 



(51) 



6e m^c3 V^ 7 W (1 + a g mjy( ^" B<,) )(l - 

in terms of and _E 7 . eo denotes the electric constant, a s and a* singlet and triplet scattering 
length, r s and r t singlet and triplet effective range, fi p and fi n the magnetic moment of proton 
and neutron, respectively. 

With this dependence of cjbl on the incident photon energy E 7 , we derive a^d^np) as the 
mean cross section of photo-disintegration per photon with E 1 > Bd'. 
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and (E,y) is given by eq. fH3|) . leading to 

_ 3 

r(n P ^d 7 ) _ 5.755 7] ( kT \ 2 b± a ( 7e ^ e7 ) 

^(■yd^np) 1 + ^ \m N c 2 ) cr^ d ^ np) 

Next, we analyse the term on the right hand side of eq. f|48j) : 

r. 



\ 2 

\ Tip 

n d J n n 



(52) 



expansion 
(yd^np) 



2 t 



n d n* c cr^d^np) 



1.040 



Gg*(T N )h 3 kT e 



Ed. 



@(ye— >ey) Tip Tl n 



{Bd + kTf af yd ^ np) n n n d 



(53) 
(54) 



Both analysed terms depend on T exponentially, but at T = Tjv, when term (1521) reaches 
unity, term fl54l is still of order 10~ 2 and therefore negligible. In other words, the expansion 
rate of the universe is still dominated by the rate of deuteron photo-disintegration, when 
the primordial nucleosynthesis starts. Of course, the expansion of the universe causes the 
opening of the deuterium bottleneck, but its major influence is the reddening of the radiation, 
which enables the deuteron production to win over photo-disintegration. This simplifies our 
task drastically and we obtain the following equation for Tjv, respectively tjy'- 



5.755 r] ( kTj 



- N 



and finally the corresponding neutron to proton ratio: 



cr^d^np) \nd J n n 



(55) 



n 



n,. 



-(Ti 



N 



exp 



X n {Tf) 

where X n is given by eq. fl33|) 



90fr 3 c 5 



(56) 



The neutrinos have decoupled from equilibrium before the annihilation of electron 
positron pairs. Therefore the entropy due to this annihilation is transferred exclusively 
to the photons, i.e. g*(T N ) = 2 + 



I 2 JV„ 



li- 



Assuming neutron conservation after tjf, we finally calculate Yp, the primordial 4 He 
abundance by weight. Since 4 He is not further transformed into heavier nuclei, because 
elements with nucleon mass number A=5 and A=8 are insufficiently stable to function suc- 
cessfully as intermediate products for nucleosynthesis at the available densities, we derive: 



^n n m He 



\n n m He + (n v 



Tin )m T . 



1 + 2 



(57) 
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where m# e denotes the mass of the helium nucleus and — (Tjy) is given by eq. fl56l) . 



This analytic expression for Yp reproduces the observation based results and the nu- 
merical results (see section results for details), assuming neutron conservation when 52 — 66 
% of the neutrons have reached stable states inside light nuclei. Therefore eq. (|57j) provides 
our basis for finding the dependence of Yp and the possible deviation of v from its present 
value t>o, m order to finally constrain ^ by recent results on the primordial 4 He abundance. 



Expressing all key-parameters of Yp by terms of v, we start with the most important 
one, the deuteron binding energy. Within our narrow range of interest (l 2 ^ 2 -) < 0.5%), 
we determine its dependence on v by yarying the pion mass in di fferent nucleon-nucleon- 
potential-model calculations, based on f Arenhovel fc Sanzone! 1991 ). and derive: 

B d {v) ~ B d {v ) (A — (A — l)v^7^), (58) 

where A is a model dependent constant as follows: 

Bonn-A potential fMachleidt et al. 1987) : A = 2.3, 

Paris potential fLacombe et al. 1980) : A = 28.2, 

Argonne V14 potential fWiringa et al. 1984) : A = 61.4. 

For our purpose, we take the mean average A = 30.6. 

As Bd changes, E 1 and the cross sections o"( 7 d_ mp ) and <J( n p-*dri) change, accordingly. 
Furthermore, we have to consider, that the mass of the electron varies proportionally 



m e (v) = m e (v ) — , 
v 



(59) 



which enters the Klein-Nishina cross section. 



Concerning r n , the mean lifetim e of the free neutron, we use the expression (IGaBner fc Lesch 



2008 


), based on ( 


Mtiller et al. 


2004) 



r n {v) ~ r n (v ) (1 - 4.88 



V . 



V 



(60) 



Next, we have to consider the change on Q, the neutron to proton mass difference, 
which influences the freeze -out concentration. We separate the electromagnetic contribution 
(IGasser fc Leutwylerlll982l ) and obtain 

v 



Q 



-0.76 + 2.0533317 



[MeV]. 



(61) 



The Fermi coupling constant Gp is related to v by (IDixit fc Sherlll988l ): 



Gp(v) 



>V2 



GeV 2 



(62) 
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Finally, we derive a relation between Yp and v, to constrain the permitted variation of 
the Higgs vacuum expectation value by the primordial 4 He abundance: 



AY P 



-38 



Av 



2.08 

v 



{N u -3\ f&G\ [An 

+ 0.0355 ( — v — — J + 0.0874 ( — J + 0.0042 In ( — - 

Varying each parameter separately (assuming the others fixed), we derive: 



/ AB d \ ( Ar n 

0.106 I I + 0.056 ' 



V B 



d 



( AGp . 
0.235 {-q 1 - 1 - : > 7)2 



/Am e \ ( Aq A \ f AN U \ 

0.006 I ] - 0.195 ( — — I + 0.0355 ( I + 



+ 0.0874 



m e 

AG 
~G~ 



V 9 a J 



V 3 



0.2602 ( ^jpj + 0.0042 In + 1 



(63) 



AQ 
Q 



(64) 



We briefly note, that using eq. ( 1631) . one can also obtain constraints on N u , G, h and rj. 



Results 



Constraining the possible time variation of v , we use ob servation based result s as well 
as the standard big bang nucleosynthesis code, developed by IWagonerl (119731 ) and iKawano 
(119921 ) . This standard code still se ems adequate for our purp ose, although newer nuclear 
reaction rates have been evaluated (jDescouvemont et al.ll2004l ). 

The baryon to photon ratio is given by r] = (273.9 ± 0.3) 10~ 10 Q b h 2 ( iSteigmanl 120061 ). 



where Qb is the present ratio of the baryon mass density to the critical density and h is the 
present value of the Hubble parameter in units of 100 kms~ 1 Mpc~ 1 . We take 100 Qbh 2 = 



2.273 ± 0.062, the 5-year mean value of WMAP (IKomatsu et al. 



20081 ). and obtain: 



6.049 • 10~ 10 < r] < 6.403 • 10~ 10 . 



Implementing rj and r n = 878.5 s (jSerebrov et al.ll2005l ). the numerical code delivers 



Y 1 = 0.2467 ±0.0003 

and thus eq. (163]) constrains the possible time variation of v. 

Av 

< 1.5 ■ 10" 
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Using the observation based results of llzotov et al.l (120071 ) we derive 

Av 



Y 2 = 0.2516 ±0.0011 



< 5.6- 10" 



and the results of iPeimbert et al.l (120071 ) lead to 



K = 0.2477 ± 0.0029 



Av 



< 1.4 -lO -3 . 



Combining the observations of H II regions (Y 2 and Y 3 ) and the numerical simulation based 
on WMAP results (Yi), according to 



AY = Yp(H II regions) - Yp(WMAP + SBBN), 

we alternatively derive (using eq. (163]) ) two more conservative estimates: 

Av 



(65) 



Y 2 -Y 1 
Y 3 -Y 1 



0.0049 ±0.0011 



0.0010 ±0.0029 



v 

Av 



v 



(2.4 ±0.6) • 10~ 3 
(0.5 ± 1.4) • 10~ 3 . 



We avoid the term "observational results" because the cited publications more or less 
consist of interpretations of the observational 4 He abundance plus theoretical input and 
constraints by the cosmic microwave background. The different interpretation as a result of 
the badly understood systematics lead to incompatible data. For consistency, w e only cite 
data based on recent He I recombination coefficients by iPorter et al.l (120051 . 120071 ). 



4. Conclusions 

Big-bang nucleosynthesis offers the deepest reliable probe of the early universe. Its 
predictions of the light element abundances play a major role in constraining cosmological 
models. The increasing precision of observational results on primordial abundances opens 
new scientific fields that can be tested by BBN. We present a calculation, how observations on 
primordial 4 He may provide insight into the fundamental property of elementary particles: 
the Higgs vacuum expectation value. We find constraints on its allowed time variation 



Av I 



< 1.5 ■ 10" 4 . 



We thank D. Tekle and A. Bauer for help on the Kawano code and S. Winitzki and A. 
Jessner for valuable comments on the manuscript. 
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